Accepted to ApJL: June 19, 2012 

Preprint typeset using I^T^X style cmulatcapj v. 5/2/11 



o 

(N 
C 

o 

(N 



6 

CO 

o3 



> 
Os 

o 

(N 



IS 



EVIDENCE FOR ASYMMETRIC DISTRIBUTION OF CIRCUMSTELLAR MATERIAL 

AROUND TYPE IA SUPERNOVAE. 

Francisco Forster 1 , Santiago Gonzalez-Gaitan 1 , Joseph Anderson 1 , Sebastian Marchi 1 , Claudia P. 
Gutierrez 1 , Mario Hamuy 1 , Giuliano Pignata 2 , Regis Cartier 1 

1 Departamento de Astronomi'a, Universidad de Chile, Casilla 36-D, Santiago, Chile and 
2 Departamento Ciencias Ffsicas, Universidad Andres Bcllo, Av. Republica 252, Santiago, Chile 
Accepted to ApJL: June 19, 2012 

ABSTRACT 

We study the properties of low-velocity material in the line of sight towards nearby Type la Supernovae 
(SNe la) that have measured late phase nebular velocity shifts (i> ne b)j thought to be an environment- 
independent observable. We have found that the distribution of equivalent widths of narrow blended 
Na I Dl & D2 and Ca II H & K absorption lines differs significantly between those SNe la with 
negative and positive v no b, with generally stronger absorption for SNe la with v nc h > 0. A similar 
result had been found previously for the distribution of colors of SNe la, which was interpreted as 
a dependence of the temperature of the ejecta with viewing angle. Our work suggests that: 1) a 
significant part of these differences in color should be attributed to extinction, 2) this extinction is 
caused by an asymmetric distribution of circumstellar material (CSM) and 3) the CSM absorption is 
generally stronger on the side of the ejecta opposite to where the ignition occurs. 

Since it is difficult to explain 3) via any known physical processes that occur before explosion, we 
argue that the asymmetry of the CSM is originated after explosion by a stronger ionizing flux on 
the side of the ejecta where ignition occurs, probably due to a stronger shock breakout and/or more 
exposed radioactive material on one side of the ejecta. This result has important implications for both 
progenitor and explosion models. 

Subject headings: supernovae: general — distance scale 



1. INTRODUCTION 

Type la Supernovae (SNe la) are important tools for 
understanding the evolution of the Universe because 
of their high luminosities and light-curve homogeneity, 
which led to their st andardization f or cosmological dis- 
tance measurements (jPhillipsI 11991: IHamuv et all fl99"?il) 
and the discovery of the acceleration o f the Universe 
(jRiess et al.lll998HPerlmutter et al.lll999l) . They are also 
important because of their complex nucleosynthetic out- 
put and high ejecta kinetic energies (10 51 erg), which 
make them key ingredients for galaxy evolution theory. 

Unfortunately, we still lack a clear understanding of 
the nature of SN la progenitors. Although it is generally 
accepted that their progenitors are carbon oxygen white 
dwarfs (CO WD) in mass transfering binary systems (for 
a review, see iHillebrandt &: Niem cvcr 2000j), there is no 
agreement on what their companions are, whether the 
CO WD reaches the Chandrasekhar mass at ignition, 
what the mechanism that triggers the ignition is and how 
the CO WD burns to form the ejecta. 

Perhaps the two scenarios that best match observed 
S N la spectra and li ght curves are those considered 
in iRopke et al.l (|2012t ): either a stably accreting car- 
bon oxygen white dwarf (CO WD) that reaches the 
Chandrasekhar mass (Mch), undergoes a thermonu- 
clear runaway and burns in a deflagration to detonation 
transition, in the Mph-single degene rate (SD) scenario 
(|Nomoto et al.lll98i IKhokhlovl[l99l : or a pair of CO 
WDs that merge and ignite dyna mically in the viole nt 
double degenerate (DD) scenario (jPakmor et alJ[2uTl . 

Two significant differences between both scenarios are 
the typical central densities at ignition and the distribu- 



tion of circumstellar material (CSM) around their pro- 
genitor systems. In the Mch _ SD scenario the central 
density at ignition is sufficiently high to avoid the syn- 
thesis of radioactive 56 Ni in favor of stable iron group ele- 
ments (IGEs) in the innermost regions of the ejecta. The 
Mch _ SD scenario has potentially abundant CSM that 
could be produced by either the wind of the donor star, 
weak nova explosions experienced by the accretor, or by 
mass loss due to inefficient mass transfer. In the violent 
DD-scenario, conversely, the primary central density is 
never high enough to avoid the synthesis of radioactive 
56 Ni in the innermost regions of the ejecta, and since 
the merger and ignition happen in a dynamical time- 
scale there is no time to leave a significant imprint in the 
CSM. 

The Mph-SD scen ario appears to over-produce IGEs 
(|Nomoto et alj 11984) , which are not produced signifi- 
cantly in the violent DD scenario. Observations of flat- 
topped profiles of Fe II lines in the optical and NIR at 
late phases has been argued to be evide nce for stable 
IGEs in the innermost regi ons of the ejecta (|Hoflich et al.l 
I2004t iGerardv et ail 120071 ). Moreover, significant nebu- 
lar velocity shifts (v ne b) measured in these flat-topped 
profiles have been suggested as evidence for off-center 
ignition sce narios, which is expe cted in Mch _ SD scenar- 
ios (see e.g. iNonaka et aLl [2012). The relation between 
u nc b and other properties such as the rate of evolution 
of the velocity of Si absorption lines (the velocity gradi- 
ent) or the colors during the photospheric phase of evolu- 
tion tiMaeda et al.ll2010l:lMaeda et al.ll201ll; iCartier et ail 
1201 ll) seems to lend support to this picture. 

For the nearby SN 201 lfe, there is evidence for lack of 
interaction between the ejecta and a companion star or a 
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companion star wind. The radius of the co mpanion star 
has b e en constrained to be l ess than 0.1 dBloom et al.l 
[20121: INueent et all l20TTt IBrown et al.l 120111), and its 
mass loss, less than 10~ 8 yr -1 (|Horesh et al.l 
120121) . However. [Sternberg et all (|20lll) showed that SNe 
la show an excess of blue-shifted narrow Na lines, which 
is evidence for CSM around type la progenitors ejected 
before explosion. Pa rticular cases, specifically SN 20 02ic 
(IHamuy et all 120031) . SN 2006X (IPatat et al] [20071 c.f. 
Chugail 120081)" SN 1 999cl (IBlondin et al.l I2009T) and SN 
20071e ([Simon et al.l l20(39h have also provided evidence 
for CSM around SNe la, although in the cas e of SN 2002ic 
the n ature of the explosion is unclear (jBenetti et al.l 
l200fft. 

ISimon et al.l (|2009f ) noted that SNe la with variable 
narrow absorption lines t e nd to h ave broad lines and high 
velocity gradients. iFolev (2012a) recently suggested that 
higher velocity at peak, redder SNe la are preferentially 
found in lower mass host galaxies, but also th at they have 
an exc ess of blue-shifted absorption systems (|Folev et al.1 
l2012b[ ). These observations highlight the importance of 
understanding the relation between progenitor systems 
and dust properties for cosmological distance determina- 
tions fe.g. lGoob ar 2008: Folatelli et aLllMtCft . Since SNe 
la with f n eb > are known to have higher velocity gra- 
dients, generally associated to higher velocities at peak , 
and to be redder (jMaeda et al.ll2010HMaeda et al.ll201lD . 
we investigate the relation between the presence or ab- 
sence of narrow Na I Dl & D2 and Ca II H & K absorp- 
tion lines (hereafter Na and Ca lines), a proxy for ma- 
terial in the line of sight, with v nc h, a geo metrical proxy 
that should depend only on viewing angle (jMaeda et al.1 
120101 ) and that should therefore be independent of any 
evolutionary or host galaxy effects. 



automatically defined. The position of the nodes is cho- 
sen to avoid being too close to the lines being measured, 
and the flux associated to every node is defined by locally 
averaging the observed flux with typical velocity ranges 
proportional to the separation between nodes and a con- 
tinuous smoothing function. After the nodes are placed 
and the continuum defined, we visually inspect the spec- 
tra, divide the observed flux by the continuum and com- 
pute the area below unity within 300 km s _1 of every 
line in A: 



EW = dX 



dX, 



(1) 



where fi are the fluxes of pixels with velocities within 300 
km s _1 of the line of interest, fi is the continuum value 
for every pixel, and dX is the wavelength separation be- 
tween pixels. In order to obtain more robust EW deter- 
minations, for every SN we use the median EW among 
all available epochs, excluding those epochs with signifi- 
cantly larger errors. 

We have considered four sources of error in our method 
somewhat conservatively: flux measurement uncertain- 
ties (<7f), continuum flux uncertainties (cr c ), possible sys- 
tematic deviations in our method (c sys ), and the median 
absolute deviation (a robust error measurement) of the 
EWs for different epochs for every SN (ctmad)- c/ and a c 
were obtained using the quoted flux uncertainties for ev- 
ery spectra at different wavelengths and propagating er- 
rors using equation (TTJ to compute the EW uncertainty, 
ctew : 



°ew — 



fi 



fi fi 



dX J 



(2) 



2. SAMPLE SELECTION AND DATA ANALYSIS 

In this work we have used the sample of 25 nearby 
SNe la with measured n ebula r velo city shifts (t^ s h) de- 
rived by IBlondin et all (|2012t ) and iMaeda et ail (|2011l ) 
for which we have measured equivalent widths (EWs) of 
Na (24 SNe) and Ca (23 SNe) absorption systems. The 
EWs were obtained from mid-resolution spectra in the 
public Center for A strophysics (CfA) spectra dat abase 
(IBlondin et al.ll2012t)Jha et al.lll999UBranch et alJ l2003t 
iMatheson et aJJl200aTWang et al.l l 2009D for most SNe la, 
from iSilverman et al.l pOHl) for SN2009dc. from high- 
resolution spectra in the European Southern Observa- 
tory (ESO) archive for SN 2000cx SN 2001el and SN 
2006X, from iSternberg etall (l201ll) for S N 2007sr and 
SN 2007af, using lYaron fc Gal- Yarn! (|2012t ). values from 
the literature for SN 1986G, SN 2004eo and SN 2006dd 



(iPhillipsl 119931: iPastorello et all 120071: IStritzinger et al.1 
[2010D, and our own X-Shooter ( jVernet et al.ll201lD data 
for SN 2010ev. We note that although there are public 
spectra for SN 1990N and SN 2004eo, we do not use them 
because they appear to have been artificially smoothed 
in the wavelength direction, removing all evidence for 
narrow lines. 

2.1. Equivalent width determination 

We have written a Python script to measure EWs in ev- 
ery SN for every epoch available. The script first defines 
the continuum by tracing straight lines between nodes 



where <jf i are the flux uncertainties in every pixel, aj i 
is the uncertainty associated to the continuum at every 
pixel, defined as the average node flux uncertainty of the 
closest two nodes to every pixel, and where each node 
uncertainty is defined by weighting the Of i 's in the same 
way as the flux of every node is defined. 

To determine cr sys we have run our script assuming 
that the lines are in positions of the spectra where no 
low velocity absorption lines should be found, changing 
these positions artificially along the spectra and using 
quadratically averaged values of the derived EWs as the 
final errors. To determine ctmad we simply computed 
the median absolute deviation of the measured EWs at 
different epochs, excluding those values with relatively 
large errors. 

2.2. Validation of the method 

We have compared our measurements of the EWs with 
those in the literature and with those obtained from pub- 
licly available high-resolution spectra, as shown in Fig- 
ure [TJ Our results are consistent within errors in 14 out 
of 15 measurements, with the only exception being SN 
1998bu, for which we obtained a significantly higher EW, 
that upon visual inspection appears to be due to con- 
tinuum definition uncertainties. We have also compared 
our measurements with those obtained with IRAF's splot 
routine and found consistent values in 90% of the cases, 
without any highly significant exception. The nebular 
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EW(Na) [A] 

Fig. 1. — Comparison of measured equivalent widths of blended 
Na absorption features due to intervening material at the distance 
of the SNe la. We show the values obtained with our Python 
script using m id-resolution s p ectra , EW(Na), and compare them 
with those in iBlondin ct al. (2009(), circles, Taubcribcrgc r~et al.l 
(|201U) . square, and high— resolution spectra measurements, stars. 
Red filled symbols are SNe la with v nch > and blue open symbols, 
SNe la with u nc t, < 0. The dotted line shows the EW(Na) comp = 
EW(Na) relation. The only SN la with significantly different mea- 
surements in our sample is SN 199 8bu, with i> ne b < 0, and for 
which we measure a higher EW than Blondin ct al. (2009). 



TABLE 1 

SUPERNOVAE WITH AVAILABLE NEBULAR VELOCITIES AND THEIR 
MEASURED EQUIVALENT WIDTHS (EWs) FOR NA AND CA LINES IN 
A AS DESCRIBED IN THE TEXT. NA INDICATES THAT DATA IS NOT 
AVAILABLE. 





Manic 


v nch [km s x ] 


EW(Na) 


EW(Ca) 


SN 


1986G 


-639 


± 600 


3.60 


± 


0.01 


1 


75 


± 


0.01 


SN 


1990N 


-478 


± 600 






NA 








NA 


SN 


1994D 


-2220 


± 220 


0.05 


± 


0.22 





06 


± 


0.25 


SN 


1994ae 


-1204 


± 600 


0.36 


± 


0.16 





33 


± 


0.22 


SN 


1995D 


-1193 


± 600 


0.32 


± 


0.21 





23 


± 


0.29 


SN 


1997bp 


2539 


± 410 


1.81 


± 


0.21 


1 


12 


± 


0.26 


SN 


1998aq 


-1157 


± 129 


0.19 


± 


0.16 





35 


± 


0.20 


SN 


1998bu 


-1351 


± 305 


0.63 


± 


0.12 





07 


± 


0.16 


SN 


2000cx 


-172 


± 600 


-0.07 


± 


0.08 


-0 


01 


± 


0.01 


SN 


2001el 


980 


± 206 


0.66 


± 


0.04 





67 


± 


0.10 


SN 


2002bo 


2258 


± 361 


2.48 


± 


0.18 


1 


20 


± 


0.26 


SN 


2002dj 


2021 


± 414 


1.09 


± 


0.33 





78 


± 


0.31 


SN 


2002er 


582 


± 600 


1.77 


± 


0.17 


1 


62 


± 


0.31 


SN 


2003du 


-443 


± 283 


0.07 


± 


0.21 


-0 


02 


± 


0.22 


SN 


2003hv 


-2825 


± 384 


-0.09 


± 


0.25 


-0 


11 


± 


0.41 


SN 


2003kf 


1070 


± 201 


0.41 


± 


0.20 





69 


± 


0.25 


SN 


2004dt 


-2623 


± 600 


0.23 


± 


0.23 


-0 


02 


± 


0.36 


SN 


2004eo 


-170 ± 25 


0.00 


± 


0.20 








NA 


SN 


2005cf 


-123 ± 600 


0.04 


± 


0.21 





23 


± 


0.20 


SN 


2006X 


2832 ± 343 


2.06 


± 


0.04 





80 


± 


0.06 


SN 


2006dd 


-1553 ± 201 


3.22 


± 


0.08 


1 


00 


± 


0.10 


SN 


2007af 


831 ± 350 


0.39 


± 


0.02 





28 


± 


0.17 


SN 


2007sr 


1185 ± 600 


0.29 


± 


0.01 





09 


± 


0.26 


SN 


2009dc 


-1190 ± 600 


0.94 


± 


0.15 


2 


00 


± 


0.31 


SN 


2010ev 


2145 ± 225 


0.67 


± 


0.05 





37 


± 


0.02 



velocities and equivalent widths used in this analysis are 
shown in Table 1. 

3. RESULTS 

In order to test whether SNe la with negative and posi- 
tive Wncb have different distributions of intervening mate- 



rial in the line of sight, we have performed Kolmogorov- 
Smirnov (KS) tests of the cumulative distribution func- 
tion (CDF) of the measured equivalent widths of Na and 
Ca lines caused by intervening material in the Milky Way 
and in the SN host galaxies. We have also tested whether 
SNe la with negative or positive v nc h have different ex- 
tinction or colo r properties using E(B-V) values obtained 



with SNooPy jlBurns et al 



from SiFTO (|Conlev et al 



2011 and (B-V) Bmax values 
" 2008ft best-fitting models. 
The KS test results are shown in Table 2. 

As expected, the distribution of intervening Milky Way 
material between SNe la with negative and positive v nc h 
is consistent with being drawn from the same parent pop- 
ulation. However, the likelihood that the distribution of 
intervening material in their host galaxies is the same is 
significantly smaller. Given that the separation of the 
sample by the sign of v ne h should be purely geometrical, 
we should not expect any significant differences in the 
distributions of host Na and Ca lines, unless the inter- 
vening material is of circumstellar nature and with an 
asymmetric distribution. 

To test whether the different distributions of EWs are 
caused by an unlikely combination of host galaxy prop- 
erties, we compare the distributions of host galaxy mor- 
phologi cal types and inclin ations from the Asiago cat- 
alogue (|Barbon et al.l Il999f) and the higher resol ution 
pixel based statistic NCR (|James &: A nderson 2006}) that 
quantifies any association with regions of star formation 
(SF) using GALEX NUV (23 out of 25 SNe) and our 
own Ha maps (14 out of 25 SNe), separating between 
SNe la with v ne b > and v ne b < 0. We do not find 
any significant difference between these two populations 
for global (morphological type and inclination) or local 
(NCR statistics) host galaxy properties. The KS test 
p- values are shown in Table 3. 
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Fig. 2. — Equivalent widths of Ca vs Na lines due to intervening 
material at the redshift of the SN host galaxies. We show SNe la 
with u no b > as red filled symbols, and those with i> ne b < as 
blue open symbols. SNe la occurring in pixels with NCRntjv > 
are shown as stars, those with NCRnuv = are shown as circles 
and those with no measured NCRnuv are shown as squares. The 
three SNe la with v nch < at high EWs are SN 1986G (open blue 
square), SN 2006dd (open blue star) and SN 2009dc (open blue 
circle). SN 2004eo, for which we only have EW(Na) information, 
is shown below the dotted line. 
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TABLE 2 

KS TEST FOR NA AND Ca LINES. THE NULL HYPOTHESIS IS THAT SNE IA WITH NEGATIVE AND POSITIVE V ncb HAVE THE SAME DISTRIBUTION 
OF COLORS OR EWS FOR A GIVEN LINE DUE TO THE MlLKY WAY (MW) OR INTERVENING MATERIAL AT THE DISTANCE OF THE SNE. 



KS test 


E(B-V) SNooPy 


( B - v )i»r 


MW Na 


Host Na 


MW Ca 


Host Ca 


p- value 


0.006 


0.002 


0.268 


0.013 


0.402 


0.030 



TABLE 3 

KS TESTS FOR ENVIRONMENTAL DIFFERENCES BETWEEN SNE IA WITH NEGATIVE AND POSITIVE ?Veb- THE NULL HYPOTHESIS IS THAT THE 
DISTRIBUTION OF GLOBAL (MORPHOLOGICAL TYPE AND INCLINATION) AND LOCAL (NCRnUV and NCRhq) HOST PROPERTIES 

DISTRIBUTIONS ARE THE SAME. 



KS test Morphological type Galaxy inclination NCRnuv NCRho 

p-value 0.133 0.390 0.768 0.944 



In Figure [2] we plot the measured EWs for the Na and 
Ca lines divided by the sign of v nc b- Red filled sym- 
bols are SNe with v nc b > and blue open symbols, 
those with v nc b < 0. Stars indicate SNe la occurring 
in pixels associated to SF (NCRnuv > 0), circles indi- 
cate SNe la occurring in pixels with no association to 
SF (NCRnuv = 0) and squares indicate SNe la with no 
NCRnuv measurements. SN 20 04eo, which only has a 
reported non-detection of Na in (|Pastorello et al.ll2007f) 
and for which we associate a typical error of 0.2 A, is 
shown below the dotted line. In this figure, SNe la with 
Vncb < appear clustered at relatively low EW values, 
with three exceptions at high EWs, whereas v nc b > SNe 
la tend to have higher EWs in general. This suggests that 
the intervening CSM absorbs more light opposite to the 
direction of where the ignition occurs. 

The fact that some SNe la show high EWs and v nc b < 
(3 out of 14) is expected, since in some cases the EWs 
will be due to host galaxy material and not CSM ab- 
sorption. Even if their high EWs are due to CSM our 
KS tests rule out that these exceptions are significant. 
In particular, SN 1986G and SN 2006dd appear to be lo- 
cated behind prominent dust lanes of their host galaxies 
and SN 2006dd has the second highest NCRnuv statis- 
tic in our sample (0.827), implying a possible associa- 
tion with SF. However, SN 2006dd and SN 2009dc have 
been suggested to show some signs of CSM based on their 
Na line evolution (IStritzinger et al.ll2010T ) and light curve 
shape (jTaubenberger et al.l 120 111 ) hence we cannot rule 
out that their high EWs are due to CSM. 

Most SNe la with u nc b > in Figure [5] have some de- 
gree of association with SF (NCRnuv > 0), which could 
be associated to higher EWs due to non-CSM in the 
line of sight. However, only 2 out of 10 SNe la with 
Vncb > are strongly associated with SF (NCRnuv > 
0.5), whereas 4 out of 15 SNe with v nc b < are strongly 
associated to SF. This suggests that higher EWs in our 
sample are not primarily due to host galaxy non-CSM, 
which is consi stent with the excess of blue-shifted Na 
lines found bv lSternberg et all (|2011| ). and which is sup- 
ported by the high p-values for the NCRnuv an d NCRh q 
KS tests in Table 3. It could also be argued that NCR 
statistics are only a measure of relative SF within galax- 
ies. To account for this, we have removed early-type 
galaxies from our sample and repeated the KS tests from 
Table 2, obtaining significantly low p-values of 0.022 and 



0.018 for the distribution of EW(Na) and EW(Ca), re- 
spectively. 

4. DISCUSSION AND CONCLUSIONS 

Our main conclusion is that the distribution of EWs 
from narrow blended Na and Ca lines differs significantly 
between SNe la with negative and positive v nc b- Because 
the sign of v ne b should be purely geometrical, it should 
not correlate with any host galaxy or average SN prop- 
erties. This suggests that part of the lines are formed by 
CSM ejected by the progenitor before explosion, which 
is found asymmetrically distributed after explosion. The 
differences in color found for SNe la with different v nc b 
( Ma eda et aLll2011[ ). which we confirm using SNooPy and 
SiFTO light curve fitting, should be explained in part by 
extinction due to an asymmetric distribution of CSM, 
and not by different ejecta temperatures with viewing 
angle alone. 

There are two possibilities to explain this result: 1) the 
CSM material is asymmetrically ejected from the system 
before explosion and is aligned with the side of the ejecta 
where ignition occurs, or 2) the CSM material is initially 
spherically symmetric and is affected by an asymmetric 
distortion after explosion, which is stronger on the side 
where ignition occurs. 

One possible physical process that could break the 
symmetry of the WD and have an effect on both the 
direction of ignition and the distribution of CSM is 
rotation. However, rotation should produce cylindri- 
cally symmetric systems, where positive or negative v ne b 
should lead to the same average CSM properties in the 
line of sight. Another possibility is the gravitational field 
of a companion star, but the t ypical velocity of f sets fo und 
for the absorbing material in lSternberg et aLl (|2011h im- 
plies that the material was ejected at a time before ex- 
plosion much longer than the expected orbital period of 
the progenitor systems, again producing some cylindrical 
symmetry in average. 

Thus, we conclude that only 2) is possible to explain 
our results. One way to produce an asymmetric distribu- 
tion of dust, EW(Na) and EW(Ca) is to have an asym- 
metric ionizing field after explosion. This would destroy 
most of the dust and ionize Na and Ca in the direction 
where ignition occurs. A comparison with p re-explosion 
imag es of the core-collapse SN 20 12 aw (|Fraser et al.l 
120121) supports the idea of dramatic changes in extinc- 
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tion after explosion. The fact that the correlations are 
found to be stronger for dust (colors) than for EW(Na) 
and EW(Ca) is suggestive, since they are in ascending 
order of ionization potential. One problem with this in- 
terpretation could be the constraints in radio and X-ray 
for CSM-ejecta interactions. This would be solved if the 
material were placed far enough from the ejecta, and only 
in some cases allow the ejecta to interact with the CSM, 
as was seen in SN 2002ic, SN 2006X and 20071e, or if the 
measured EWs are due to saturated lines in low mass, 
clumpy CSM with a wide distribution of velocities. 

Our results provide hints about the explosion asymme- 
try and CSM properties of SN la progenitors. If a larger 
sample confirms that most SNe la with v nc h > have 
high EWs of Na and Ca absorption lines, the majority 
of progenitors should contain significant CSM. If this is 
the case, most SNe la with u nc b < should have more 
strongly ionizing fields (higher energy photons) after ex- 
plosion, either during shock breakout, or afterwards if 
more radioactive material is exposed in the ignition side 
of the explosion. 

Finally, significant CSM is c onsistent with som e Mchr~ 
SD or slow DD scenarios (e.g. iFrver et al.ll2010l ) and not 
with violent DD scenarios. However, in our inter preta- 
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